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ABSTRACT 


This thesis discusses the theory behind the gas dynamic laser. 
A computer program that solves the equations is included along with 
typical outputs. A description of a test apparatus is also presented 


with proof of laser action. 
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I. INTRODUCTION 


The operation of any type of laser system depends on the pro- 
duction of a population inversion. Descriptions of the phenomena 
that take place in a laser may be found in [1] and [2]. There are 
several methods of producing a population inversion, one of which is 
by means of a gas dynamic process. The gas dynamic method usually 
incorporates a combination of rapid heating and rapid cooling of a 
gas. 

In the rapid cooling phase, a gas that is initially at a high 
temperature and pressure is exapnded by flow through a supersonic 
nozZle, which causes a sudden decrease in the temperature and pressure. 
The rapid drop in temperature is accompanied by a decrease in the 
population of the upper energy levels, each of which decays with a 
Characteristic time constant. In certain gas mixtures the time 
constants are such that a lower energy level depopulates much more 
rapidly than a higher energy level. The result is a population 
inversion of these two states. The method has been successful in 
producing a laser output as shown by Kuehn and Monson [3]. In their 
experiment a large gas reservoir was heated and pressurized. The gas 
was then passed through a nozzle into an evacuated dump tank. Mirrors 
were placed in the sides of the nozzle to form the optical cavity. 

It was the purpose of this thesis to build and test a gas dynamic 
laser utilizing a shock tube assembly to generate the initial high 
temperature conditions of reservoir gas, and to simulate the device on 


the digital computer according to the theory developed by Basov [4, 5]. 
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II. THEORETICAL DISCUSSION 


The fact that a gas dynamic laser is theoretically possible was 
established by Basov, et al. [4]. This work presented the require- 
ments that must be met by a set of energy-level relaxation times if an 
inversion is to be established. The cases of rapid heating and rapid 
cooling were discussed. The ground states of a CO, -N, mixture were 
found to satisfy the general requirements. Basov, et al. [5] examined 
this system in detail and offered a set of equations to describe the 
phenomenon. Part of the results of their work is a system of equations 


that describe vibrational relaxation. These equations are 
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This gives a population inversion of 
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In the equations, nitrogen is the "qa" gas and carbon dioxide is the 


"8" gas. The oo represent the number of collisions between a specific 


molecule of type i and any j type molecule. Poo is the probability of 


resonant collision between a co, and an No molecule. Pop is the corres- 


ponding function for two co, molecules. oF m are the probabilities of 


> 


transition from one quantum state to another. The n and m correspond 
to the different energy levels shown in Figure l. A listing of the 
(n,m) that goes with a given transition may be found in [5]. The 0. 


are the characteristic temperatures of the corresponding levels, and 
hv, 
are given by Q. = rs . Here, hv. is the energy of the given state 


above the ground state. The above relations give the rate of change 


of energy in specific energy levels of CO, and No: EB through ee 


represent the energies stored in the specific energy level. The initial 


0,/T_1) 


values of E, through Eq are set by the relation E, = (RO, /M, )/(e i 


1 
[5 5 mca 

A basic explanation of the origin of these equations starts by 
considering the change of population of some state of one particle as 


the particle undergoes collisions with another particle. The rate of 


change of population would be given by 
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where one is now considering an N-level system. Pom is the probability 
that a particle will undergo a transition from the nee energy level to 
th lea jean 
the m energy level. (Basov uses Pam ©° indicate a probability of 
transition due to resonant collisions and Qj to indicate a probability 
due to other than resonant collisions.) The first term on the right- 
: ; 2 710 
hand side gives the depopulation of the i level. The second term 
: ; sth ; 
introduces the repopulation of the i level. Ze is the number of 
> 
collisions that an e type particle will undergo with a p type particle. 
Note that only transitions due to interaction with another particle 
are being considered at this point. Transitions that may be caused 
by other reasons (spontaneous, photon stimulated, etc.) are not included. 
These were assumed negligible. 
In order to describe the population inversion process, one must 
know the transition probabilities (pee) and the collision rates Z ‘ 
? 
The actual calculations that must be carried through to obtain the Poe 
will not be worked out here. Rather, an intuitive feel for the method 
used will be presented. One starts by writing Schroedinger's equation 
for the collision process [6]. This will result in a wave-function 
description of a particle that is under the influence of an intermole- 
cular interaction potential. This interaction potential is a function 
of the distance between the involved particles. 


The resulting wave functions can be written down as some combination 


of the wave functions of an isolated particle. The set of constants 


ess 





that perform this transformation is directly related to the transition 
probabilities Dees The relationship between the elements of the set 
of constants and the transition probabilities is known, and is of 
sufficient simplicity as to allow calculation of the probabilities. 
The origin of the a "s is discussed in [5] where they are derived 
from the kinetic theory of rigid molecules. The probabilities given 
by Basov are altered versions of the theoretical relationships so that 
the best experimental correlation results. 

In order to form a closed system of equations that describes the 
gas dynamic process, one must describe the temperature, pressure and 
density as functions of time. This is necessary due to the dependence 
of the probability relationships on these parameters. The quasi one- 
dimensional flow equations form a set of independent equations that 
describe the above variables. One of the relationships is the equation 


of state: 
R 
P= o0 M a (5) 
where R is the universal gas constant and M is the molecular weight of 


the gas. If there is more than one gas, this equation can be separated 


into the partial pressure components. That is, 


Assuming steady-state flow, another relationship is the continuity 


equation, which is given by 


* (pav) = 0. (6) 


A third relationship is that for the conservation of momentum, as given 


by 
Vio dp 
PY ax dx (7) 
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Finally, there is the conservation of energy for adiabatic flow. This 
equation takes the form: 

GR +e )+8 Skreet oven (8) 

pO MY Q a B 2 0 
and in effect states that there is only a finite amount of energy that 
can be distributed among the different modes of energy storage. 

These eight equations allow one to describe the relaxation pehnomenon 
as exhibited by an adiabatically cooled gas. A digital computer program 
that solves the equations is shown in the Appendix. 

Before describing the theory of the gas dynamic process, keeping 
the theory of the N, -CO, gas dynamic laser in mind, the N, -CO, mole- 
cular system will be discussed. The 10.6-micron N,-C0, molecular laser 
system inverts the population of certain of the molecular vibrational 
energy levels of the ground state of the co, molecule. The vibrational 
energy levels are split by the molecular rotational energy. However, 
in the energy-level diagram, for clarity purposes, these rotational 
energy levels are not shown. The different vibrational modes of the 
ground state of the co, and No molecules are shown in Figure l. 

With reference to Figure 1, ER is the lowest vibrational level of 
the N, atom. Since No is diatomic, there is only one vibrational degree 
of freedom. This is a symmetrical stretching mode that can be visualized 

as similar to two balls that are attached by a spring. 

(00°1) is the energy of the assymetrical stretching mode of the CO, 
molecule. This is the case of the two oxygen atoms moving in one 
direction while the carbon atom moves in the other’direction. 


1 
(01 0) is the bending mode energy. This is a degenerate state, as 


the molecule can bend in either of two perpendicular planes, the degeneracy 
being indicated by the superscript. 
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(10°0) is the energy of the symmetric stretching mode. This is 
the same as in the case of N,, except that there is a "stationary" 
carbon atom between the vibrating oxygen atoms. 

The above energy levels are simple harmonic oscillator states, and 
therefore the energy difference between the nth and (nth+l) energy 
levels is the same, regardless of the value of n. The nomenclature 
for the higher vibrational energy levels is the same as above. For 
example, the level (02°0) is the nondegenerate bending mode that has 
an energy of approximately twice the (0110) level. Note that although 
it is not shown, the (02°0) would be right below the (10°0) level. 

There are other levels in the region shown by Figure 1, but the 
levels shown are those that enter into the strong 10.6-micron transition, 
and are the most important ones for the present purposes. 

With the energy levels set forth, it is now necessary to explain 
several other phenomenon before the laser can be completely described. 
The first of these phenomenon is what is known as a collision of 
the second kind or a resonant collision. Basically the resonant 
collision is the interaction of two atoms, one of which is in an 
excited state. The other atom may be in any state. The result of 
the interaction is that the excited atom transfers its energy to the 
second atom, while returning itself to the ground state. This may be 


written schematically as [2]: 





xor Y* Ske a Y 





where the * indicates an atom in an excited state. It has been shown 
that the probability that the interaction will occur decreases as the 
energy gap between the energy level of atom Y and the closest available 


energy level of atom X increases. If the energy difference is greater 
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than 800 er the collision will effectively not occur [7]. Although 
a proof will not be presented here, some insight can be gained by 
considering the fact that both the energy and the momentum of the two 
atoms must be conserved. This criterion is most readily fulfilled when 
there is little "internal" energy (i.e., electronic, vibrational, or 
rotational) that must be transformed into kinetic energy of trans- 
lation of the atoms. 

The next item that must be discussed is the Ny energy-level diagram. 
As seen from Figure 1, for No there is only one state of interest: the 
lowest vibrational level of the ground state. This state is a molecular 
level, and since N, has no net electric dipole moment, there can be no 
transitions directly to the ground state due to any electromagnetic 
interaction. The result is that for the Ny molecule to return to the 
ground state from the vibrational levels it must undergo a collision 
with another particle (other molecules or the container walls). This 
is known as a metastable state, and is characterized by a long lifetime. 


With the above information, the CO molecular gas dynamic laser 


a2 
system can be described. Once the gas has been heated, populating the 


upper levels, it is then expanded by means of a supersonic nozzle. 


The excited No molecules interact with the CO, molecules, and the two 


types of molecules undergo resonant collisions which serves to maintain 


the population of the 00°1 level of the CO The population of the 


2° 


other CO, energy levels depopulate rapidly as both levels (shown in 


2 
Figure 1) are connected to the ground state by large transition 
probabilities. The large energy gap between these states and EO pro- 


hibits repopulation due to resonant collisions. This creates a 


population inversion between the 00°1 and the 10°0 states. There are 


ig 





other inversions, for example between the 00°1 and 02°0 levels, but 
the above inversion is the strongest, and gives the 10.6-micron 
transition. The lasing action may now occur. 


There is one other characteristic of the N,-CO, laser system that 


should be pointed out. This is the fact that any impurities in the 
gases will cause a marked decrease in the output power. This is due 


to the resonant transfer of energy from the Ny molecules to the 


impurity atoms, which robs energy that would have otherwise been 
pumped into the inversion. 


Although the above theory does describe the same N, -CO, laser 


system, higher output powers can be obtained by the use of additives. 


In the case of a He -N, -CO, laser, a major factor in the population of 


the 00°1 level appears to be due to cascading effects from higher 
energy levels of both molecules [8, 12]. That is, there is resonant 


transfer of energy between the CO, and He at higher energy levels 


Z 


(collisions of the second kind at upper levels). The co, molecules 


: Oo : : 
then spontaneously drop into the 00 1 level to set up the inversion. 
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IIIT. EXPERIMENTAL APPARATUS 


The previous section dealt with the theory behind the gas dynamic 
laser. The apparatus used to test the theory is shown in Figures 2 
through 13. The shock-tube facility used had been previously built 
at the Naval Postgraduate School and it is discussed in detail by 
Penaranda [13]. 

Figure 2 shows a schematic of a double-diaphragm shock tube. The 
shock tube generates a shock wave by suddenly expanding a high-pressure 
gas into a low-pressure gas. The shock wave may be visualized as 
formed by the superposition of infinitessimal compression waves. The 
shock wave is accompanied by large gradients in temperature and pres- 
sure. The double-diaphragm arrangement was used because of the repro- 
ducability of run conditions that it affords. In the double-diaphragm 
assembly, one drops half of the driver pressure across each diaphragm. 
The diaphragms (flat brass plates) are cut to withstand this pressure. 
However, when valve #1 is opened the entire driver pressure is felt 
across diaphragm #1 which has been cut to rupture well below this 
pressure. Therefore, the shock will always be originated by the same 
total driver pressure. This system eliminates the necessity that the 
diaphragms rupture at an exact pressure. The only requirement is that 
they rupture in the pressure range between the full driver pressure 
and one half that value. 

Figure 3 shows the regions of interest for the case of a reflected 
shock wave. The shock wave is traveling into the driven gas (region 1) 
at a velocity S. The contact surface is propagating into the driven 


gas at a velocity C. This surface separates Region 2 from Region 3. 
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Regions 1 and 2 contain pure driven gas while Region 3 contains pure 
driver gas. That is, there is no mass flow across the contact surface. 
The array of lines in Region 4 represents the rarefaction wave, which 
is produced by the expansion of the driver gas. Region 4 is the driver 
gas at rest. Region 5 contains pure driven gas after the shock has 
reflected from the end wall of the shock tube. The conditions of the 
gas in Region 5 were used as the initial conditions for the experiment. 
Analytical descriptions of the state of the gas in any of the 
regions may be obtained by means of the conservation of mass, momentum 
and energy equations in conjunction with the adiabatic gas laws. This 
is done in (9) and (10). The results of interest to the present experi- 


ment are 


oo 5 = [2-01] [Grp w= 2 0-0) ana 


Cats Coin 2 
Ts [2 (y-1) M+ G-y)) [G-7-1) MW - 2 G-1)) 
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where y is cp/cv, a. is the speed of sound in Region 1 and M is the 
mach number of the shock wave. The conditions in Region 5 can be 
controlled by variation of M, P, and a,la,. 

The remainder of the physical apparatus is shown in Figures 4 
through 12. Figure 4 is a schematic of the entire setup. The three 
amplifiers shown serve to amplify the signals from the thin film 
gauges (T.F.G. in the diagram). The distance between T.F.G. #1 and 
#2 is 5 ft. 3 in. Therefore by recording when the shock passes each 
of the gauges, the actual shock speed can be determined. The schematic 
of the amplifiers is shown in Figure 5. 

The biasing circuit for the detector is shown in Figure 6. The 


detector is biased at about 370 microamps and gives a signal of several 
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millivolts. The detector is gold-doped germanium with a bandpass of 
Ped to 10:6 microns. The time constant of the detector is 50 
nanoseconds. 

The "dump tank" serves as a reservoir for the gas after it has 
passed through the optical cavity formed by the mirrors. As will be 
seen in a latter portion of this section, the observation time is 
on the order of 2 milliseconds. The particle velocity will be super- 
sonic after the gas expands through the nozzle. For a particle 
velocity of 4000 fps, the dump tank must be 4 feet long to avoid 
interference of the returning gas with that gas that is still expand- 
ing through the nozzle. The gas that reflects from the end wall of 
the dump tank will introduce a change in the conditions within the 
gain path. This change in the conditions marks the end of the test 
time. The dump tank can be seen in Figure 7. The fingerlike pro- 
jection from the endwall is a safety release valve, which is set to 
release when the pressure exceeds 200 psi. This is well below the 
rated pressures of all components in the system and will therefore 
prevent any overpressures from building up. 

The mirror and nozzle assembly can be seen separated in Figures 8 
and 9. A mirror mount is shown disassembled in Figure 10. A cutaway 
view is shown in Figure 11. The head of the bolt (shown in Figure 10 
with the mirror attached) is ground off so that the mirror surface is 
recessed 0.075 inch from the nozzle wall. This is to help preserve 
the mirrored surface. Original test runs were made with gold-coated 
glass mirrors (not recessed) and mirror surface deterioration was 
severe. In some places the gold was dislodged from the glass. The 


recessing of the mirrors alleviated this problem. 
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Thesmirror has a radius of curvature of 4.5 inch (2.25-inch focal 
length). The gain path with the mirrors recessed is 3.14 inches. The 
mirrors are 98% reflective at 10.6 microns. 

The nozzle is a linear nozzle with a half angle of 30°. The throat 
is 2 millimeters high. The center of the mirrors are 1.0 inch down the 
nozzle from the throat. The inside of the nozzle is roughened in order 
to minimize internal reflections of spontaneously emitted photons. 

The plexiglass shield shown in Figure 12 absorbs any stray 
radiation. The end wall is at 10° to the incident beam (which enters 
from the left of Figure 12). This will relfect 7% of the energy. The 
reflecting surface is an IRTRAN 2 crystal. The hole drilled in the 
plexiglass behind the IRTRAN is to isolate the IRTRAN from the load 
(that is from another piece of plexiglass that will absorb the trans- 
mitted beam). The two metallic bars on the far side serve as counter- 
balances to the detector (which is filled with liquid Ny under operation). 
The white wire carries the signal while the black lead grounds the 
detector and Dewar flask. 

Figure 13(a) shows the output of the thin film gauges of the shock 
tube for a test run (0.5 m sec/cm). The conditions were 600 psi Ny 
driving 40 mm Hg. ARGON was the driven gas. The resulting mach speed 
was 4.0. An endwall pressure gauge was substituted for the nozzle and 
its output is shown in Figure 13(b) (0.5 m sec/cm). The flat region in 
the middle of the trace is the observation time for Region 5, and is 


found to be on the order of 2 m sec. 
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IV. RESULTS 


A. THEORETICAL 

Figures 14 through 18 show different outputs from the computer 
program. Figure 14 shows a comparison of the program solutions to 
the predictions by Basov, et al. [5]. The case considered was for 
an axially symmetric nozzle, with initial test conditions of 10 
atmospheres and 1500° Kelvin. The two solutions were found to differ 
by no more than 3%. 

Figure 15 shows a comparison of the linear nozzle as used in the 
experiment to the hyperbolic nozzle considered by Basov. It should 
be noted that the pressure drops faster in the hyperbolic nozzle. 

The effect that this has on the population inversion can be seen in 
Figure 16. The two larger curves show comparisons of the computer 
predictions to Basov's predictions for the hyperbolic nozzle. The 
smaller population inversion in the linear nozzle is caused by the 
higher pressure and temperature (as compared to the hyperbolic nozzle). 
Also, since the temperature decreases more slowly in the linear 

nozzle, the peak in the population inversion occurs farther down the 
nozzle. 

Figure 17 shows the population inversion at the center of the 
optical cavity as a function of temperature. There is some optimum 
temperature for each initial pressure and co, concentration. Below 
this temperature less energy is pumped into the upper levels. Higher 
temperatures increase the gas velocity through the nozzle. Since the 
characteristic relaxation times are effectively independent of tempera- 


ture in this range, the population inversion peak occurs farther down 
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the nozzle than the center of the optical cavity. This causes a 


smaller inversion at the center of the optical cavity. 


B. EXPERIMENTAL 

Figure 18 shows the result of one of the initial tests. Time 
limitations precluded further experimentation. The setup used in 
this initial test for lasing action was a modification of that shown 
in Figure 4. The modified setup used only the dual-beam oscillo- 
scope. The oscilloscope was triggered by the first thin-film gauge. 
The second and third gauges were fed into different (i.e., opposite 
polarity) differential inputs of one channel. This was the lower 
trace in Figure 18. The detector output was fed directly into the 
second channel. The positive jump in the lower trace at 3.8 msec is 
the difference of the outputs of the second and third thin-film 
gauges, and comes just prior to reflection of the shock wave off the 
endwall. At this time the detector output shows a negative step, 
which was caused by the lasing action. Two milliseconds later the 
conditions of the gas were upset and the observation time as described 
by the computer program was over. This is noted by a slight ripple 
in the upper trace at 6.0 msec. The large spikes are due to recreation 
of favorable conditions by internal shock reflections. Their large 
size is due to mixing of the helium driver gas. This was shown by the 
fact that the spikes disappeared (but not the original step) when Ny 


was used as a driver gas. 
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V. CONCLUSIONS 


From the reproducibility of the output from the detector, it was 
concluded that lasing action was taking place. However, insufficient 
data were taken to permit comparison of the actual laser output to 
the theoretically predicted value. The only indication that the theory 
correlates with the experiment is the fact that a population inversion 
was predicted for the test conditions that resulted in laser output. 

Further work with the system, as shown in Figure 4, is necessary 
to supply the data required for full comparisons between the theory 
and experimental results. Further experimentation with the laser 
should supply insight into both the validity of the theory and into 
the shock reflection and gas mixing that caused the large output 


spikes. 
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APPENDIX 


The eight equations given in.Section I form the closed set 
necessary to describe the gas dynamic laser. The computer program 
at the end of this section solves these equations. Descriptions of 
the variables may be found in tables I and II. There are a few 
portions of the program that need clarification. 

The factor (TT/TC) in statement 25 ensures that T = TT during the 
fuest iteration. Lf T # TT, a divergent solution results. In some 
cases the solution may converge after several iterations. The error 
between T and TT is small and is due to round-off errors of the 
conversion factors used. 

The factor (0.98) in the relation for CSQ (c?) also ensures 
computational stability. That is, in order to compute dp/dx, one must 
divide by Gee / coe At x = 0 (at the throat), v = c, and one would 
therefore be dividing by zero. The factor 0.98 ensures that this 
situation does not occur. The resulting error is less than 1%. 

The subroutines RKLDEQ and OSPLOT may be found in the USNPGS 
System 360 program library. More information on them may be obtained 


there. 
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Variable 


¥(1),¥(2),¥(),¥(4) 
Hower, V(5) 


P1AM 


TE, El, TS ,T,1C 

Y(6), VT 

CSQ 
F(1),F(2),F(3),F(4) 
F(5) 

F(6) 


Er 


TABLE I 


Physical Interpretation 


Total energy in given level 
Pressure 


Actual driven pressure 


Temperatures 

Gas particle velocity 
Speed of sound squared 
Population rates 
Pressure rate of change 
Velocity rate of change 


Normalized population 
inversion 
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Units 


ant es 
atmospheres 


millimeters 
mercury 


degrees Kelvin 
cm/sec 
pee eae 

-1 
cm /gr/sec 
atm/sec 
cm/sec/sec 


dimensionless 





TABLE II 


Computer Variable Basov Variable 





' 
¥(1) ,Y@),¥(3),¥(4) E,> Ey", Ey, E, 
Y¥(5) P 
Y (6) Vv 
F(i) FiC1) 
TS T 
oO 
P5 P 
Oo 
CONC Y 
Z.: Gon 
it 1j 
ore 855 
Bere | 
ij 1J 
PI AN/N 
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2330.7 
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Ground State 
COs No 


FIGURE | 


CO> AND No ENERGY LEVELS (in cm!) 


Si 





Volve 


# | 


End 
Wall 


+2 # | 
Diaphragms 


FIGURE 2. 
DOUBLE DIAPHRAGM SHOCK TUBE. 
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Shock Tube At Time t, 70 
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Time 
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Region 3 





Centered 
Rarefraction 
Wave 








Regior: ~ 


Distance 


FIGURE 3 


PLOT OF DISTANCE VS. TIME IN SHOCK 
TUBE AFTER DIAPRAGM IS BROKEN. 
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Shock Tube End Wall 
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ADJUST NULL BIAS WITH A VTVM BEFORE 
‘ CONNECTING OUTPUT TO AMPLIFIER. 


3 lOOK Pot Adjusted 
For Nuli Bias 
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FIGURE 6 
BIASING . CIRCUIT 
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EIGURE =” 
END VIEW OF SHOCK TUBE 
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FIGURE 8 
MIRROR MOUNT 
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FIGURe 9 
NOZZLE ASSEMBLY 
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FIGURE 10 
MIRROR MOUNT (DISASSEMBLED) 
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FIGURE It. 
CUTAWAY OF MIRROR MOUNT. 
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FIGURE 
PLEXIGLAS SHIELD AND DETECTOR 
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OUTPUT OF T.EG. FOR TEST RUN 
(0.5 m.sec./ cm) 





FIGURE 13(b) 


OUTPUT OF ENDWALL PRESSURE GAUGE 


(0.5 m sec. /cm) 
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Detector 





2 msec/cm 


Lower and upper trace: 10 mv/em 


Run conditions: 
Driven gas: 200 mm (Hg) of 20.7% cO,, in Ny 


Driver gas: 600 psi of 40% Argon in Helium 


FIGURE 18 


DETECTOR OUTPUT SHOWING LAS#R ACTION 
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